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An effect of paramagnetic impurity located in a vicinity of a quantum well (QW) on spin polar- 
ization of the carriers in the QW is analyzed theoretically. Within approach of Bardeen's tunneling 
Hamiltonian the problem is formulated in terms of Anderson- Fano model of configuration interaction 
between a localized hole state at Mn and continuum of heavy hole states in the InGaAs-based QW. 
The hybridization between the localized state and the QW leads to resonant enhancement of inter- 
band radiative recombination. The splitting of the configuration resonances induced by splitting of 
the localized state in magnetic field results in circular polarization of light emitted from the QW. 
The developed theory is capable of explaining known experimental results and allows for calculation 
of the photoluminescence spectra and dependence of integral polarization on temperature and other 
parameters. 

PACS numbers: 75.75.-c, 78.55. Cr, 78.67.De 

I. INTRODUCTION 

Various phenomena based on interference of a discrete quantum-mechanical state and a continuum of states have 
been intensively studied since the famous paper by U. Fanoi rated among the most relevant works of 20th century^. He 
suggested a theoretical approach often regarded as Fano- Anderson model or Fano configuration interaction which was 
successfully applied to explain the characteristic asymmetric resonances observed in atomic spectroscopy experiments. 
It further appeared that numerous examples of Fano resonances existed in atomic and nuclear physics, condensed 
matter physics and optics^. The co-existence of a localized discrete level and continuum of states within the same 
energy range is also quite common in low-dimensional semiconductor structures'^—. The structures containing a 
quantum well (QW) and a ferromagnetic or paramagnetic layer located in the vicinity of the QW nowadays are 
gaining a special interest as they are believed to combine the high mobility of the carriers in the QW and magnetic 
properties provided by the magnetic layer. In particular, an exchange interaction with ferromagnetic layer leads to 
spin polarization of holes^. For GaAs-based structures with a J-layer of Mn the holes probably play an important 
role in promoting ferromagnetic state of the Mn layer^^. The system considered in the present work consists of GaAs- 
based hctcrostructure with a In^jGai-xAs QW {x = 0.1 — 0.2) and (5-laycr of paramagnetic acceptor (Mn) located at 
a distance of a few nanometers from the QW. A number of recent experiments show that the Mn (5-layer gives rise to 
circular polarization of the photoluminescence from the QW in magnetic ficld^ii^. It was found that the intensities of 
the radiation at the wavelength corresponding to interband direct transitions in the QW differ for opposite circular 
polarizations in magnetic field applied normally to the QW plane. If Mn is replaced by a non-magnetic acceptor 
(carbon) no polarization is registered. Thus the observed polarization is not due to g-factor of the free carriers in the 
QW which in this way is proved to be small. On the contrary, the holes localized at Mn do have g « 3ii and possibly 
can penetrate into the QW by means of quantum-mechanical tunneling which is expected to be of resonant-type 
if the energy of the localized state coincides with that of free 2D hole in the QW. We use the Fano configuration 
interaction approach to study tunneling hybridization between a localized hole state at paramagnetic impurity with 
the 2D continuum of states in the QW. We show how this hybridization should reveal itself in the photoluminescence 
at the QW wavelength. The theory allows for calculation of circular polarization degree in magnetic field. We have 
to mention here that there are other mechanisms that might contribute to the observed polarization like non-resonant 
tunneling of electrons from QW to Mn and recombination with the holes localized at Mn. While it still remains unclear 
which of the mechanisms has larger contribution to the polarization in experiments on photolumenescence^ii^, in our 
paper we restrict our consideration by configuration interaction for the holes quite naturally expected in the p-type 
system under study. 

II. HOLE STATES AT ACCEPTOR AND IN QUANTUM WELL 

We assume that the eigenstates of the holes in the QW form a continuum of states starting from the energy of size 
quantization (taken for zero of energy scale). These 2D states can be characterized by two quantum numbers, say 
and ky are the projections of 2D wavevector on the axis x and y in the plane of the QW. Considering interaction with 
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a point-like impurity it is more convenient use cylindrical basis. In this case each state in the QW is characterized by 
the magnitude of the wavevector k and index I denoting the cylindrical harmonic. The corresponding wavefunction 
is given by: 



7^1/4 ^ 



(1) 



where Ji{kp) is Bessel function of order I, p and are the polar coordinates in the plane of the QW, S is the area of 
the QW plane, rj [z) is the envelope function of size quantization in z-direction. The wavefunction is normalized to 
unity. Firstly, let us treat this continuum as a set of discrete states characterized by energies £k and wavcfunctions 
Lpki. Considering the free 2D carriers implies the magnetic field applied to the QW to be nonquantizing. The validity 
of this assumption is discussed in section IVII Below we consider Ina;Gai_a;As QW having only one level of size 
quantization for the heavy holes, we neglect the light holes being split off due to size quantization. Thus the basis of 
Bloch amplitudes to be used is formed of the states with certain projection of total angular momentum J = 3/2 on z 
axis which is perpendicular to the QW plane: 



(e3/2, ei/2 , e_i/2 , 6^3/2) . 
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The wavefunctions V'fcij in this basis have the form: 
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The kinetic energy of a state in the QW is related to the wavevector value as 
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where is the in-plane heavy hole mass in the QW. 

In order to determine wavefunction 1/; of a hole localized at an acceptor one should consider the kinetic part of the 
Luttinger Hamiltonian and attractive potential of the acceptor U (r). The spherically symmetrical potential preserves 
the symmetry Fg. thus the ground state is 4- fold degenerate and can be classified by projection of angular momentum. 
The eigenfunctions of Luttinger Hamiltonian with spherically symmetric attractive potential can be explicitly found 
in the model of zero radius potentiaU^. In the basis of Bloch amplitudes they are expressed as follows: 
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Eq is the binding energy of a hole at the acceptor, Yim are the spherical harmonics. mih,mhh - respectively are 
the bulk masses of light hole and heavy hole in GaAs. Note, that the radial part of all nonzero components of the 
wavefunctions (jl]) have two characteristic decay lengths, the largest of the two being always determined by the light 
hole mass m;^ . 

III. TUNNELING BETWEEN ACCEPTOR AND QUANTUM WELL FOR COMPLEX BAND 

STRUCTURE 

For the studied system the potential barrier separating the Mn (5-laycr from the QW is weakly transparent for 
tunneling so the localized hole state mix up with the QW states. Rigorous calculation of eigenfunctions of the 
resulting state is rather hard to perform as it requires solving stationary Schrodinger equation in the complicated 3D 
potential and with account for the complex valence band structure. In order to circumvent the explicit solving of 
Schrodinger equation for tunneling problems the so-called tunneling or transfer Hamiltonian formalism is commonly 
used as originally proposed by Bardeeni^. Within this approach the total tunneling Hamiltonian is expressed as 
H = Hmu + Hqw + Ht, where Hmu is partial Hamiltonian having the localized hole state at Mn as its eigen 
state. Hqw in the same way corresponds to the QW itself, its eigenfunctions ^pk form non-degenerate continuum of 
states characterized by quantum number k. The term Ht accounts for the tunneling. In the secondary quantization 
representation the total Hamiltonian can be written as follows; 



where a"*", a - the creation and annihilation operators for the localized state characterized by its energy eg, and 
c^,Cfc - the creation and annihilation operators for the continuum state having energy e^. For the convenience both 
energies Eq and £k here and below are measured from the level of size quantization of holes in the QW so that 
is simply their kinetic energy. The expression ^ is rather general, in fact it can be regarded as introduction of the 
coupling between two systems into the Hamiltonian in the most simple phenomenological way. From this viewpoint 
the coupling parameter Tk is still to be determined through exact solving of the eigenvalue problem for the whole 
system. Bardeen's approach suggests a simple recipe for calculation of the tunneling parameter Tj. (also referred as 
tunneling matrix element) for the case of weak tunneling through potential barrier: 



where K is the kinetic energy operator, commonly equal to 




For this problem the 3D integration ([7|) straightforwardly reduces to integration along any surface inside the barrier. 
In order to apply the same approach to the problem of the holes tunneling in GaAs it has to be generalized for the 
case of complex band structure. It would be tempting to do it by treating K in ([7]) as the kinetic part of the effective 
Hamiltonian and ■i/'j fk as envelopes in the basis of Bloch amplitudes. For the valence band of Ina;Gai_a;As - based 
heterostructure the holes ground state with total angular momentum J = 3/2 is described by Luttinger Hamiltonian 
{Tikx ^^ky ^fik^ are, as usual, the momentum operators along the appropriate axis): 
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The functions (pk in (O become 4-component vector functions in the basis Analogously to the simple band 
case integration over the whole space is reduced to the integration over the surface 0,$ inside the barrier, at that, 
only z-projection of the kinetic energy operator is required. For the eigenfunctions in the QW with zero light hole 
components the expression for tunneling parameter simplifies into: 

n = iB-A)jjs[^,*^^^-^±^k*). (10) 

Regrettably, this straightforward generalization of ([7]) fails to be correct. Indeed, the decay length of wavefunctions 
(j4|) is determined by light hole mass while the decay length of the of the QW states ([3]) is governed by the heavy hole 
mass. Due to this circumstance the result of the surface integration ((TU)) becomes crucially dependent on the particular 
position of the integration surface inside the barrier. However, it can be shown that in the case of two masses the 
exponential dependence of the tunneling parameter Tk on the barrier thickness is determined by the smallest mass. 
The value of is also governed by a formfactor corresponding to the overlap of the wavefunctions along the QW 
plane. We consider the kinetic energy of a hole in the QW being small compared to the QW depth, i.e. k « q. 
In this case the overlap of the localized wavefunction occur only with the zeroth harmonic of ([1]). The only nonzero 
tunneling matrix elements are therefore the following: 

Tk = Tik+T2k 

Ti, = {A-B) J^^ dS (RoYooJo {kp) £rj (z)) , 
^{A-B) J^^ dS (^i?2F2o Jo (M £v (z)) , 

The calculation of these integrals presents no difficulty for any particular surface inside the barrier, however further 
detalization of these values is beyond the accuracy of the approach. The only statement valid is that the tunneling 
parameter Tk has an exponential dependence on the barrier thickness with the light hole mass entering the exponent 
index: 

n = (A - B) C {k/q) ^'^y cxp (-X (k/q) qd^) , (11) 

where 1<X^2, C^l are weak dimensionless functions of k/q. Our further analysis will be focused on the case of 

the hole's kinetic energy being substantially less then the binding energy Eq, i.e. k « q. For this case it is reasonable 

to assume that Tk does not depend on k. However, its rapidly decreasing behavior for k ^ q has to be kept in mind as 

it provides convergence for any integration over k involving Tk- For the case k « q the particular shape of the QW 

while affecting the particular value of the tunnelimg parameter does not play an important qualitative role and can 

be assumed rectangular. The interaction with only zeroth cylindrical harmonic means that the continuum spectrum 

modified by tunneling is non-degenerate. While this fact is not principal for the qualitative results obtained below, it 

simplifies the mathematics. From the analysis given above we conclude that the tunneling configuration interaction 

exists only between ipkQ _3 and s-nd the same interaction (governed by the parameter Tk (|lip ) is between ipkQ +3 

and 'ipj_a . 
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IV. THE SPECTRUM MODIFIED BY TUNNELING 



The transfer Hamiltonian ([5]) with known tunneling parameter (jlip allows one to construct the eigenfunctions of 
the whole system given those of the localized state ip and the QW states (pk- 

^{E)^VQ{E)^ + Y,^k{E)ipk, (12) 

k 

E denotes the energy of the state Here fk are the wavefunctions with zeroth cylindrical harmonic, as was shown 
above the other harmonics are not affected by the tunneling configuration interaction. Plugging (|12p into the stationary 
Schrodinger equation: 

iJ* = E^ 

with H being the effective Hamiltonian ([5]) one gets the following system of linear equations: 

I'oeo + I] VkT^. = Evc, 

k (13) 
VkEk + Tfci^o = Evk 
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Solving the eigncnvalue problem for system (jl3p one can get the spectrum and the coefficients j/q, i^fc, i.e. the eigen- 
functions of the system. Transition from discrete set of states Vk{E) to continuous function i>{e,E) is straightforward 
(as the continuum states are non-degenerate we can use the energy e instead of k as the quantum number of the 
state). Instead of ([T^ and we write: 

/•CO 

^{E)^vo{E)^+ v{E,e)ip{e)de, (14) 
Jo 



(E) £o + /o°° t (e) f {E, e) de Evq {E) , , . 

v{E,e)e + t{e)vQ{E) = Ev{E,e) . ^ > 

The normalizations for if) and ^p{e) are: 

(£o)|^(eo)> = 1, 

{^{eme'))^5{e-s'). (16) 
With the chosen normalization the discrete tunneling parameter and the one entering (|15p are related as follows: 

TVo (£) = e is) , (17) 

where 



is the density of states for the selected basis of states with zeroth cylindrical harmonic. The discrete system is an 
eigenvalue problem, but the continuous problem ^T5\i is not. In the present work we consider the case of the localized 
energy level lying within the range of the continuum: Eq >> t^. For this case the solution is obtained as shown in^: 

° nH^E) + {E-eof' 

V {E, e) = (E) (P^^ + Z{E)t {E)6 {E - e)^ , (19) 

where 

E-eo-F{E) 
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Jo [E - e) 

P stands for the principal value and Eq is the center of configuration resonance, which appears to be slightly shifted 
from Eq. 

eo{E)=eo + F{E). (21) 

Because of fc << g it is reasonable to treat t =const everywhere, except for (PD)) where decrease of t at i? oo is 
necessary for convergence of the integral. In order to analyze the influence of the configuration interaction on the 
luminescence spectra wc have to calculate matrix element of operator M describing interband radiative transitions 
between the hybridized hole wavefunction ^'(i?) and wavefunction of an electron in the quantum well of the con- 
ductance band £,kjc, here fee is the magnitude of the electron wavevector, 4 is the number of cylindrical harmonic 
analogously to ([T|). We assume that (a) there are no radiative transitions between the localized hole wavefunction tp 
and the 2D electron wavefunction £,k^i^ thus the matrix element for transitions from the localized state: 



6e 



M 



^)=0, (22) 



(b) the interband radiative transitions between the free 2D states in the QW are direct, the matrix element given by: 

Mo = (a.i. \m\ ipu) = UkS {k - k,) 6ij^ , (23) 
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where Uk is the appropriate dipole matrix element. With use of ([14)) . p9|) . ([22)) and ()23p we arrive to the the matrix 
element for transitions between states ^(-E) and ^k^o (according to previous notes this matrix element differs from 
Mo only for the zeroth cylindrical harmonic) : 



M 



'if {E)) = V {E,aee)u{aee) ■ 



(24) 



a = me/mjj^, where nif. is the effective in-plane electron mass, Se = h^kl/2me. The particular form of M ([24)) 
prevents from calculation of the ratio M'^/Mq as done in the classical Fano resonance calculationsi. The latter 
assumes unperturbed matrix element Mq to be constant which is obviously not the case for the direct transitions 
demanding the momentum conservation ([23|) . In our case the ratio M'^ /Mq doesn't readily give physically meaningful 
result due to the delta function in ([T^. one rather have to proceed to calculation of an observable. With the Fermi's 
Golden Rule for the transition probability we write: 
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(25) 



where Eg is the bandgap and hw - energy of radiated photon, fe,fh - the energy distribution functions for the 
electrons and holes respectively. To deal properly with the delta function entering in ([^^ and emerging in ([^5)) 
we pass on to averaging the W{hjjj) over a small spectral interval of width Q, centered at ujq: 
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Using JMj) and (HD) we obtain: 
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where / (£") = /e (a^^i?') fk {E'). The first term in brackets describes the transition rate for radiative recombination 
in the QW with no account for the tunneling, therefore N{E') here is the total density of states (including not only 
the zeroth but all cylindrical harmonics): 



(27) 



Integration assuming the functions eo-,t,u^f being constant within the range of integration {t{E) = t, u{E) — u are 
assumed constant everywhere) yields : 



W{hLOo)='^u^J{E) 
n 
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where 
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We assume weak tunneling, being the smallest energy scale. In the vicinity of resonance 
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expansion of ([28| to the first order in t"^ gives: 
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(31) 
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FIG. 1: (Color online) Modification of the luminescence spectrum by tunneling configuration interaction: numerical calculation 
for the discrete levels (solid line) and analytical formula (|32|l (dashed line) 

Note that (pij) has a term —\/h£l which does not depend on the tunneling. Its appearance is due to peculiarity of 
the mathematics of the Fano model reflected in (|19p . When a non- interacting state with energy eg is appended to the 
system so that eg lies within its spectrum, one of the energy levels of the whole system becomes doubly degenerate. 
This fact is not properly accounted for in ([T9l) and one state is lost. It should be added back manually to the 
spectral density by canceling the second term in ()3ip . Treating the same issue in a different way, one should examine 
/SW = W — Wq instead of W itself, Wq being the unperturbed transition rate: ([^5]) evaluated for t = 0. In a similar 
way studying the ratio of matrix elements in original Fano worki circumvents the disappearance of one level. 

The results obtained for a single impurity can be also applied to an ensemble of impurities provided their interaction 
between each other is weak compared to the tunnel coupling with the QW. If the concentration of the impurities is 
low enough to produce only weak perturbation of the luminescence spectra, we can simply multiply the tunneling 
term by the number of impurities. After normalization by the area of the QW we finally get the spectral density of 
the luminescence intensity: 

where n is the 2D concentration of impurities. The last term in brackets corrects the lost level issue to provide exact 
canceling of the perturbation of the spectra at t = 0. For high concentration of the impurities the formula p2p may 
give a meaningless result (the intensity may become negative at some points) . Indeed for high concentration the real 
physical picture becomes slightly different - interaction between the impurities splits their energy levels forming a 
small range of discrete levels, accordingly, the configuration resonances become slightly shifted. Taking this effect 
into account eliminates the puzzling behavior of (|32p at high concentration but docs not affect the answer for the 
calculation of polarization given in the next section. 

The analytical result p2p was verified by numerical simulation performed for the discrete system (jl3p . The system 
was solved for 500 discrete levels with interlevel separation 10"^ eV, the discrete tunneling parameter was taken 
Tfc = 3.3 • 10^^ eV, which corresponds to the continuous tunneling parameter = 10~^ eV, the other relevant 
parameters were: w = 5 ■ 10^^ eV, n = 10^° cm^^, nie = 0.03 mo, mhh = 0.5 toq, to'^,^ = 0.15 ttiq. In both 
calculations all the states were assumed fully occupied, i.e. the energy distribution function was kept f{E) = 1. 
Analogously to ([25)1 the matrix element for the discrete system was taken: Mk (e,ee) — UkSk.k^Su^. The calculation 
result presented in Fig [T] demonstrates perfect agreement with the analytical expression ([5^ and confirms the validity 
of the latter. The considered tunneling configuration interaction thus gives rise to the luminescence intensity within 
a certain spectral range ((30)) corresponding to the resonance. This increase is compensated by the decrease of the 
intensity outside of the this range as can be seen in FigHJ The width of the resonance is determined by f2 which has 
the meaning of spectral resolution of the measurement setup. However, for comparison with experimental spectra the 
inhomogeneous broadening should be accounted for as it usually exceeds the instrumental spectral resolution. An 
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expression for the integral intensity over the whole spectra reads: 



This formula follows from (|25p in the same way as ([^5]) and ([5T|) were obtained. The delta-function here is added 
manually to treat the lost level issue - it provides canceling of the second term in the limit < — > and thus gives the 
correct expression in the absence of tunneling: 

oo 


Note that the spectral width of the resonance ft does not enter the expression for integral intensity . 
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V. POLARIZATION OF THE SPECTRA 



While the redistribution of the spectral density does not change the integral intensity it gives rise to the integral 
polarization of the spectra in the magnetic field as illustrated by Figl^J The 2D holes with projection of total angular 
momentum j = +3/2 and j = —3/2 recombine emitting respectively right- (cr^) and left- (cr~) circularly polarized 
light. In section Hill it was shown that the heavy holes with j = —3/2 interact basically with the eigenfunction 
ijj_3 of the localized state. Let us denote the corresponding energy of this state Eq . The 2D holes with j = +3/2 
interact in turn with 4'+^ which corresponds to the energy e'^ . An external magnetic field applied in z would cause 

Zeeman splitting between Eq and Eq . The splitting A — e'^ — Eq may also originate from exchange interaction of 
holes with spin-polarized Mn ions. The value of A in this case is determined by exchange constant and depend on 
the degree of Mn spin polarization. The splitting of the localized energy level leads, in turn, to the splitting of the 
configuration resonance. Indeed, as follows from (p6|) . ((28)) . ([20]) the difference in the positions of the resonances E+ 
and E" corresponding to the localized levels e^j" and Eq is given by: 



A = E^ 



E- 




(35) 



Unless the position of the resonance E is too close to the valence band edge the last term in ([55)) can be neglected 
and A = A = — £q . The applicability of this result is limited to the case Eq > A. This condition, in fact, simply 
means that the splitting of the localized level does not bring any of the sublevels beyond the energy range of the 2D 
continuum so that the exploited Fano approach remains valid. Our consideration will be always limited to this case. 
With account of the energy distribution functions for the holes and electrons the shifted positions of the resonances 
lead to the difference in the luminescence intensity for the opposite circular polarizations. Let /"*", /" be the integral 



luminescence intensities of circular polarizations cr^ 
polarization is given by: 



and a respectively. Assuming |/* — /o| << /q the integral 



With use of dMj) this yields: 
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The slow varying functions in the upper integrals may be assumed as constants taken at Eq ,£q , the tunneling 
parameter will be treated as a constant in the whole range of interest {E) = . 
Then expanding over gives for the first order term: 
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FIG. 2: (Color online) Mechanism of polarization of the luminescence. The localized hole levels split in magnetic field. Each of 
them effectively couples with the 2D holes having certain projection of angular momentum. Shifted positions of the resonances 
with account for temperature distribution of the holes cause the difference in intensities of circular polarizations cr^, 



The formula ([57]) leaves not much room for further simplification for a general case of /e (E) , (E) being two Fermi 
distributions characterized by chemical potentials fi^^ fJ-h and the temperatures Tg and respectively, all four pa- 
rameters being different. Let us analyze a few particular cases leading to compact analytical expressions for P. All 
the cases imply Sq > A. Firstly, let the holes be fully degenerate and both energies £q and fo^ lying well beyond 
the quasi Fermi level of the holes so that: /x/j — >> A. In this case the distribution function of the holes can be 
considered as constant Fh{E) = 1 in the range E £ (e,^,ejj') Assuming the electrons to be non-degenerate with their 
temperature Te the formula (j37|) reduces to: 

p ^ P e-T^ sinh-^, (38) 
2kT* ^ ' 



where 

Pi 



Here T* = aT^. Exactly the same expression is valid for the case when both electrons and holes are non-degenerate. 
The only difference from the previously considered case is that now the effective temperature T* is given by 

The expression (|55)) is plotted in Fig |V] for different values of the parameter 7 = A/eq. The polarization shows 
nonmonotonous behavior with increasing the temperature. In the discussed theory the polarization arises from 
splitting of the configuration resonances positions for (t+ and spectra. The configuration resonance itself causes 
the redistribution of the transitions rate in the vicinity of the resonance energy conserving the total rate, thus the net 
polarization is subect to the difference in occupation numbers for ej^ and Eq^. The maximum integral polarization 
is therefore naturally expected when the derivative of the combined distribution function f{E) reaches its maximum 
value within the range E S (ej^,?^). For the considered case the maximum of the derivative is at eg when Eq = kT* 
and the value of the derivative decreases with increase of Eq. This explains the overall decrease of the maximum 
polarization with decrease of 7 in Fig|Vl For another case we consider the electrons distribution function /e being 
nearly constant within the configuration resonances. This can be due to their non-equilibrium distribution with a 
high quasi Fermi level or the electrons temperature Tg being much higher than that of the holes. The holes are now 
considered to have Fermi distribution function with the quasu Fermi level fih and the temperature T. We also assume 
kT << Eq. In this case from ([57)) we get: 

p_p f 2exp(/30sinh(C/2)+7 \ 

° Uxp(2/?e)+2exp(/30cosh(e/2) + iy ' ^ ^ 
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FIG. 4: (Color online) Temperature dependence of polarization for the case of electrons distribution function being constant 
within the configuration resonances. The parameter /3 = denotes deviation of the holes Fermi level from the configuration 

resonance, 7 = 0.1. 



where 

o _ £0 - 

^ ^ A ' 

^0 = (40) 

The dependence ((39)) of P/Pq on 1/^ is plotted in Fig |V] for different values of the parameter /3 (the value of 7 was 
taken 0.1). In this ease the maximum of the distribution funetion derivative is at the holes Fermi level /i/i, therefore 
the largest integral polarization corresponds to /3 = 0. For this particular case (|39p simplifies into: 
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P = Po (tanh(0 



2 cosh^ (0 



(41) 



The integral polarization obtained within given approach can be alternatively expressed through an effective g-factor 
of the holes (?e//- Let us consider Zeeman term in the Hamiltonian of the 2D holes: 

where is the angular momentum projection operator, /ig is the Bohr magneton and B is the magnetic field applied 
along z- direction. The polarization of 2D holes due to the splitting between the two subbands with j — +3/2 and 
j = — 3/2 is given by: 



Pb 



J{f{E^ Ab/2) -f{E + Ab/2)) dE 


JifiE- Ab/2) + fiE + Ab/2)) dE 





(42) 



where = 3^o5e//^- For the nondegenerate case one gets 



Pb 



tanh ■ 



2kT* 



(43) 



Comparing (|43l) with (p8| gives: 

9eff 



2kT* 
1 



-tanh 



2TmhH'^ \ _ H, Ao 
e fcT^ smh ■ 



mhh'SokT* 



2kT* 



(44) 



In the same way an expression for the degenerate case can be easily obtained. 



VI. DISCUSSION 



The key advantage of the Fano approach utilized in the present study is that the unknown eigenfunctions of the 
complex system are expressed through the known ones of the uncoupled states, in our case these are the hole localized 
at Mn and the free 2D hole in the QW. Given the expansion ((T^ any effects on the localized state can be easily 
translated into effect for the whole coupled system. For the polarization effect under study the key parameters are the 
energy of the localized level Eq, the splitting parameter A and the tunneling parameter t'^. The binding energy for a 
hole at a single Mn in GaAs is known to be i^o ~ 110 meV-i^. For the enhanced Mn concentrations in the delta- layer 
up to 10^'^ cm^^ the impurity band is established with the binding energy lowering down to 50 meV or even lessi^^— . 
Thus we consider the QW depth for the holes starting from 50 meV to be likely for the resonance tunneling effects 
under study to occur. Estimations for the splitting energy A subject to both exchange interaction between the hole 
and Mn and external magnetic field. For small concentration of Mn the splitting between the states having opposite 
projections of the angular momentum can be simply estimated as A = fJ-ogB with B being an external magnetic 
field and g sa 3 is the g-factor for the hole at Mn. This makes A '--^ 0.1 meV for i3 ~ 1 T. Samples with higher Mn 
concentrations up to a few percent are known to exhibit ferromagnetic properties'^, in this case the levels splitting A 
is to be considered with account for p-d exchange interactioni^. The particular value of A for typical experimental 
samples still remains questionable, in the ferromagnetic regime the splitting is probably believed to be in the range 
A ~ 1 — 10 meV by the order of magnitude. The magnitude of the tunnel coupling is, of course, the key parameter 
determining the polarization. From (|lip and (|18p it follows that the tunneling parameter can be estimated as 

t^ = {A-B)^ wi_hh_Ug (-45-) 

Substituting band parameters for GaAs B — A k, \ ■ h^/2mo, Eq = 100 meV, d = A nm one obtains the characteristic 
value ~ 0.1 meV. We then take the Mn concentration n ~ 10^^ cm^^, ^ = 1, eo ^ 3 meV, /3 = 0. This 
set of parameters gives Pq ~ 0.2, Pi « 0.6. The experimental temperature dependence of polarization obtained 
ii>iS qualitatively agrees with (|4T|) . Beside the analytical expressions for the general case ((37)) and particular cases 
P8 p .(|4i p . a numerical simulation of the luminescence spectra can be performed based on (|32p . An example of such 
calculation is shown in Fig lVIl For the calculation the following parameters were taken: A = 1 meV, n = 10^^ cm~^, 
T = Te = 20 K, £0 = /i/t = 1 meV, the inhomogeneous broadening of the spectra was accounted for by normal 
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FIG. 5: (Color online) An example of calculated luminescence spectra for right (cr+) and left {a ) circular polarizations. The 
parameters used in calculations are given in the text. 

distribution of Eg with dispersion taken as cr = 3 meV (corresponds to the fluctuation of the QW width by half a 
monolayer). The calculated spectra presented in Fig lVIl seem to be in good agreement with the experimental results 
obtained ir^i^. 

As mentioned in section |lT] the magnetic field in the developed theory was assumed nonquantizing. Indeed an 
estimate for the energy of Landau levels separation gives: 

hbJc — —, — ~ 0.3 meV 

for the magnetic field B ~ 0.5 T. This value is substantially less than the typical kinetic energy of the holes estimated 
as e sa 1 — 10 mcV. However, this value is comparable with the tunneling parameter t^. Therefore for the experimental 
data the validity of the developed theory is well justified for i3 < 0.5 T. 

VII. SUMMARY 

We have presented a theory treating tunneling configuration interaction between a continuum of states in the QW 
and a paramagnetic impurity located outside of the QW. The coupling between the localized state and the QW 
is described by means of tunneling Hamiltonian. We utilized the well known Fano approach for calculation of the 
matrix elements for direct radiative transitions between electrons and holes in the QW. At that the new results were 
obtained not covered by the conventional Fano formulaii^, the tunneling interaction between the localized hole state 
and continuum of states in the QW results in the symmetrical redistribution of the luminescence spectral density in the 
vicinity of the resonance. While this redistribution does not affect integral luminescence intensity, it causes an integral 
circular polarization of the light emitted from the QW provided the localized state is split in the projection of the hole 
angular momentum under external magnetic field or due to exchange interaction with Mn ions. The presented theory 
expresses the cigcn states of the system with weak tunnel coupling through the wavcfunctions of the hole localized 
at paramagnetic center and the 2D states of the continuum. For this reason it seems to be capable of describing 
other effects expected in such systems like an anisotropy of the holes g-factor in the QW induced by a paramagnetic 
impurity or the indirect exchange interaction between the localized hole states provided by the continuum located at 
a tunneling distance. 
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